Nonlinear Hydrodynamics from Flow of Retarded Green's Function 
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We study the radial flow of retarded Green's function of energy-momentum tensor and i?-current 
of dual gauge theory in presence of generic higher derivative terms in bulk Lagrangian. These are 
first order non-linear Riccati equations. We solve these flow equations analytically and obtain sec- 
ond order transport coefficients of boundary plasma. This way of computing transport coefficients 
has an advantage over usual Kubo approach. The non-linear equation turns out to be a linear 
first order equation when we study the Green's function perturbatively in momentum. We consider 
several examples including Weyl^ term and generic four derivative terms in bulk. We also study 
the fiow equations for _R-charged black holes and obtain exact expressions for second order trans- 
port coefficients for dual plasma in presence of arbitrary chemical potentials. Finally we obtain 
higher derivative corrections to second order transport coefficients of boundary theory dual to five 
dimensional gauge supergravity. 
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I. INTRODUCTION AND DISCUSSION 

Fluid/gravity correspondence has become an interest- 
ing aspect of current theoretical physics research after 
beginning of RHIC program in 2000. The experimental 
data implies that the quark-gluon plasma (QGP) pro- 
duced in RHIC are in a new state called thermalized 
matter. The evolution of QGP and hadronic matter in 
this state can be described by hydrodynamics. The tem- 
perature of the gas of quarks and gluons produced at 
RHIC is approximately 170MeV which is very close to 
the confinement temperature of QCD. At this high tem- 
perature they are not in the weakly coupled regime of 
QCD. This new phase of nuclear matter is known as the 
the strongly coupled quark-gluon plasma (sQGP). Ob- 
viously to study their properties the usual perturbation 
theory does not work. One needs different approach to 
deal with these strongly coupled system. 

The holographic hydrodynamics (or fluid/gravity cor- 
respondence) is an important tool for understanding 
some properties of strongly coupled CFTs in terms of 
the dual AdS black holes physic£] [I] -[25]. 

Low frequency (long wavelength) fluctuations of any 
interacting quantum field theory should be described by 
hydrodynamics. Low energy behavior of strongly coupled 
field theory (with gravity dual) is governed by a weakly 
coupled black hole space-time in one higher dimension. 
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QCD is approximately conformal at sufficiently large energies. 



On the other hand in classical theory of gravity a black 
hole can be viewed as a fictitious fluid membrane (with 
hydrodynamic characteristic) living on the horizon of the 
black hole |15] . Therefore from UV/IR point of view it 
would be interesting to understand the precise relation 
between the membrane fluid and the low frequency de- 
scription of strongly coupled field theory sitting at the 
boundary of AdS black hole space-time. 

One can read off the transport coefficients like shear- 
viscosity coefficients, of boundary plasma from its re- 
tarded Green's function of stress tensor 



^xy,xy\"'lJ-) 



I / dtdxe'''-'^{[T,y{x),T,yiO)]) . {1.1) 



We write the Green's function in powers of momentum 
and it is given bjrj 

(1-2) 
where, p is spatial dimension and p > 3 and 

rj : shear viscosity coefficient, 
Ttt : relaxation time for shear viscous stress, 
K : a new coefficient defined in Ref[9, 10] 

and id is frequency and q is spatial momentum (k = 
{co,0,0,q}). 

In |27| it was shown that, at low momentum (fc^ — > 0), 
the evolution of the boundary Green's function is inde- 
pendent of the radial direction and hence, it can be com- 
puted either at horizon or at boundary. Computing the 
Green's function at black hole horizon one can show that 
the shear viscosity coefficient, which is a first-order trans- 
port coefficient of the boundary plasma, can be obtained 
from the low frequency characteristics of the membrane 
fiuid. In [S5H5DJ , this issue has been generalized for grav- 
ity theory with generic higher-derivative interactions in 
arbitrary backgrounds. It has been proved that the mem- 
brane fiuid does give the correct shear-viscosity coeffi- 
cient of the boundary plasma in arbitrary higher deriva- 
tive gravity and just the knowledge of near-horizon geom- 
etry of the dual AdS black hole is enough for computing 
the shear- viscosity coefficient [5T| . 

To specify the boundary plasma completely one also 
needs to understand its higher-order transport coeffi- 
cients. For this one needs to move away from the low 
frequency (fc^ — >• 0) limit. In this case, the Green's func- 
tions flow non-trivially with the radial direction and the 
flow depends on full black hole geometry. Although the 
boundary plasma and the membrane fiuid have same 
shear-viscosity coefficients, other transport coefficients 
can certainly differ and it is not clear how the two are 
related. In Fig. [TJwe plot the radial evolution of response 
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FIG. 1. Flow of Green's function from horizon to boundary 
for two derivative gravity. 



function for two derivative gravit}]^ In zero frequency 
limit flow of the response function x is trivial. Its real 
part is constant (the constant is one in our scaling) and 
imaginary part is zero. For finite uj and q the response 
function (both real and imaginary part) has a non-trivial 
evolution. But the horizon value of the function for non- 
zero frequency is same as the horizon value of the function 
for zero frequency. Therefore we conclude that for two 
derivative gravity dual, the full momentum response at 
the horizon automatically corresponds to only the zero 
momentum limit of the boundary response. 

However in presence of higher derivative terms in the 
action the situation is different. The full momentum re- 
sponse at horizon depends on spatial momentuirr] We 
plot the flow of response function for Weyl'^ interaction 
in Fig. [2] 
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FIG. 2. Flow of Green's function from Horizon to boundary 
for higher derivative gravity for different values of lo and q. 

From this plot we see that in presence of higher- 
derivative interaction the horizon value of the response 
function depends of fc^ unlike two derivative theory. 

We also see that the imaginary part of x diverges at 
boundary. This is usual UY divergence. One has to add 
a counter term to cancel this divergence. 

In this paper, we study the flow equations for retarded 
Green's function of boundary theory analytically and 



Here we have dropped the frequency independent part of the 
Green's function. 






^ We define the response function x — 

* We will also see this analytically in sections III 
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find higher order transport coefficients of the boundary 
plasma solving this equation. We generalize the analy- 
sis for generic higher derivative gravity theory and also 
for R charge black holes. The Sow equation for Green's 
function is a first order non-linear differential equation of 
Riccati type. Because of its non-linear nature it is hard 
to solve this equation exactly. After a change of vari- 
able one can reduce this non-linear equation to a second 
order linear homogeneous differential equation. But to 
solve this we need to specify two boundary conditions. 
In this paper we deal with the non-linear equation and 
specify the boundary condition at the horizon. There- 
fore the hydrodynamic characteristic of the field theory 
at UV fix point is determined by IR boundary condition. 

For two derivative Einstein-Hilbert action the flow 
equation of retarded Green's function has been derived 
in [27]. But it is not obvious how to generalize the flow 
equation for higher derivative gravity. The derivation 
given in |27j was based on the canonical form of gravi- 
ton's action. In this paper we have considered generic 
higher derivatives terms in the bulk Lagrangian. Follow- 
ing the prescription given in j29) we construct an effective 
action for transverse graviton which has the canonical 
form in presence of any higher derivative terms in the 
bulk, and derive the flow equation for Green's function. 
Solving this flow equation perturbatively in uj and q we 
obtain second order transport coefficients namely t^ (re- 
laxation time) and k of the dual plasmar] In this way 
of computing the transport coefficients has an advantage 
over usual Kubo approach. In Kubo approach, one has to 
first find the transverse graviton by solving a second order 
differential equation and then compute regarded Green's 
function. Instead, the ffow equation is a first order dif- 
ferential equation (although non-linear). As we want a 
perturbative expansion of Green's function in powers of 
UJ and q the equation turns out to be a linear first order 
differential equation. Thus, technically, it is simpler to 
get results for causal hydrodynamics, particularly when 
the dual bulk theory is complicated. 

The paper is organized as follows. We have worked 
in five-dimensional bulk theory, the dual gauge theory 
is four-dimensional. In section |Tl] we review the deriva- 
tion of flow equation of boundary Green's function for 
two derivative gravity. Then solving this flow equation 
we compute second order transport coefficients k and r^ 
of boundary plasma. Our results matches with previ- 
ous computations of [3 [TU]. In section III we present 



In section 

the fiow equation in presence of generic higher derivative 
interaction. We calculate the higher derivative correc- 
tion for T^ and k in section |IV[ We concentrated mainly 
on Weyl^ and four derivative correction. While our re- 
sults for Weyl'^ correction are in agreement with results 



of [35], the four-derivative corrections are new results of 
this paper. We also compute the effect of gauss-Bonnet 
to T^ and K. From our results we obtain the bound on 
Gauss-Bonnet coupling constant. In section |V] we con- 
sider R charge black hole in bulk and find the second 
order transport coefficients in presence of finite chemi- 
cal potential. We also find the higher derivative effects 
on transport coefficients of field theory plasma dual to 
gauge supergravity theory. Finally, we analyze the flow 
equation of Green's function for boundary R current in 
section |VT] Appendix \K\ and [B] compliments some discus- 
sion on boundary terms. In appendix [C] |D] and [EJ we 
have provided some long expressions. 



II. FLOW OF RETARDED GREEN'S 

FUNCTION OF ENERGY-MOMENTUM 

TENSOR 

In this section we briefly review the work of Liu and 
Iqbal ^27j. They considered leading Einstein-Hilbert (E- 
H) action with a negative cosmological constant in 4+1 
dimensions and studied the motion of a transverse gravi- 
ton in this backgrouncFJ The action is. 
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g{R+l2). 



The background has a black-brane solution as, 
dS — gttdt + grrdr -\- gijdx dx-' , 



(II.l) 



1 



9tt 



5r 



4r2(l - r2) 



(11.2) 



The solution is asymptotically AdS and it has a boundary 
topology R X R^. The horizon of the space-time is at 
r — 7> 1 and asymptotic boundary is at r — > 0. 

We study the graviton's fluctuation in this background. 



9xy = g^^y^ + Kyir, x) = gl^J [1 + e$(r, x 



(o)r 



(11.3) 



By plugging it in the action and keeping terms to or- 
der e^, we obtain the following effective action for the 
perturbation 



S 



1 



167rGs 



^ ' p,q=0 



(II.4) 

Here we use Fourier transform to work in the momentum 
space k = {— w, fc}. 



<i>(r, x) = 



d'^k 



'^{r,k) 



(11.5) 



^ From Weyl invariance one can show that there are other trans- 
port coefficients in second order hydrodynamics 1101 . However 
from the expansion of retarded Green's function it is only possi- 
ble to compute only two of them. 



^ Here we restrict ourselves to five space-time dimensions, but the 
discussions are quite generic and can be extended to arbitrary 
dimensions. 



and (j)^P^ (r, k) denotes the p*^ derivative of the field 
(j){r, k) with respect to r (p + q < 2). 

Next, wc integrate by parts to obtain the bulk action 
in the following form (up to some total derivative terms) 



S 



cf^k 



dr{Ai{r,k)(l)'{r,k)(f,'{r,~k) 
+A(r,fc)^(r,fc)^(r,-fc)), (11.6) 



where. 



Ai{r,k) = - 



1 rr 

2y ^ 
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IGttGs 



(II.7) 



From this action, we can find the conjugate momentum 
n(r, fcp) of the transverse graviton (for r-foliation) and 
the equation of motion, 



and 



U{r,k^) = 2Ai{r,k)cl)\r,k) (Il.i 



n'(r, k^) - 2 Ao{r, k)cl){r, k) = . (11.9) 



The on-shell action reduces to the following surface 
ternj^ 



s-E 



d'^k 



{Ai{r,k)^'{r,k)cf>{r,-k)). (11.10) 



Following the AdS/CFT prescription given in [5], the 
boundary retarded Green's function is given as, 



Gnik^) = lim 



2A(r,fc)0'(r,fc)0(r,-fc) 
(j)o{k)(j)f){-k) 



(11.11) 



where, 4>o{kfi) is the value of the graviton fluctuation at 
boundary. Full solution of the graviton can be written as 
(j){r, fc^) = (po{k^)F{r, fc^), where F{r, k^) goes to identity 
at the boundary. We can rewrite the boundary retarded 
Green's function as, 



Gnikf,) = lini 



n(r,fc^) 



r^o 4>{r, kfj,) 



Let us define a response function of the boundary theory 



X{k,^,r) 



n(r,fc^) 

iuj(f>{r, k^) 



(11.13) 



"^ We will discuss about other boundary terms in appendix PbI 
* We set the zero frequency part of G to zero, as it gives contact 
terms 



where w = k^. This function is defined for all r and k^. 
Therefor the boundary Green's function is given by. 



Gnik^) = lim iujx{k^,r). 



(11.14) 



We will study the radial evolution of the response func- 
tion x{kfi) from horizon to boundary. Differentiating 
equation ( |II.13 ) and using the equations of motion (II. 9) 
we get, 



drXik,i,r) = iu 



9rr 

9tt 



Xikf,, 



T(r) 



S(r) 



OJ^ 



where we define 



E(r)--2A(r,fcp 
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T{r) = 2Ao{r,k^)J- 



9tt 



9r 



(11.15) 

(11.16) 
(11.17) 



Putting values of Ai and Aq given in (|IL7|) we can 
easily recover the fiow equation given in [2 7 PI However 
for future requirements, here we present it directly in 
terms of the coefficients of the graviton action. 



As mentioned earlier, the flow equation in (11.15) is 



valid for any value of momentum. This is a first order 
differential equation and we need to specify one boundary 
condition to solve this equation. That naturally comes 
from the behavior of the equation at the horizon. De- 
manding the solution to be regular at the horizon, we 
get the following condition. 



Xik,„r) 



E(r)T(r) 



w^ 



(11.18) 



For two derivative gravity this boundary condition im- 
plies thalpl 



X{k^,l)^- 



S(1)T(1) 



1 
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(11.19) 



which is independent of k^. Therefore the full momen- 
tum response at the horizon corresponds to only to the 
zero momentum limit of boundary response, x(fcp, 1) = 



(11.12) x{k^.^O,r~^0)^ 



With this boundary condition, one can integrate out 
the differential equation ( II.15[ ) from horizon to asymp- 
totic boundary and obtain the AdS/CFT response for all 



^ just notice that gtt there is negative of what we used here. 
^^ Wc choose the negative brunch. The sign of the boundary condi- 

tion lll.lSl depends on the choice of coordinate. In our coordinate 

the boundary is at r — > hence we need to choose the negative 

branch. 
^^ However, in higher derivative gravity we will see that the 

x(fc(j, 1) depends on spatial momentum. 



momentum k^. In particular, it is trivial to see that at 
(oj, ki) — >■ limit, the flow is trivial 



drX{kp.,r) = 



(11.20) 



and using the boundary condition (11.18) we get the first 
order transport coefficient of boundary fluid, i.e. the 
shear viscosity coefficient coefficients turns out to be 

In this paper, we will go away from (w, fcj — >■ 0) limit. 
As we have already mentioned, it is possible to integrate 
the flow equation for any momentum (perturbatively) 
and we can easily flnd the higher order transport co- 
efficients. The usual Kubo approach to compute these 
coefficients requires the full profile of the transverse gravi- 
ton in black hole background background (solving a sec- 
ond order differential equation), where as, using the flow 
equation, one can get these transport coefficients without 
explicit knowledge of the graviton's profile. 



A. A renormalized response function 



When we solve the flow equation (11.15 1 to get the 



boundary response function in general it involves diver- 
gence at the boundary (r -^ 0). These are usual UV 
divergences and to remove them we need to re-normalize 
the response function properly. 

We follow the holographic renormalization prescription 
of [3H[35]. As the graviton is massless, we only need to 
add the following counterterm to the graviton's action. 



Sn = 



16nG 
In momentum space, 
1 r d^k 



^-— f d''xy^^^e,x)a^e,x). (11.21) 



5, 



c 



j(l}{d,k){g'*iV'+kik')4>{6,~k). 



64^G5 Jr^s (27r)4 

(11.22) 
Therefore the renormalized Green's function is given by. 



Gr = lim 



n(?',fc^) y/^ 

0(r, k.) 32ttG5 



{g'^Lu^ + k,k') 



(11.23) 



However we will study the flow of un-renormalized re- 
sponse function defined in (11.13 1 and we define our renor- 



malized response function as, 

X^^"(r,fc,,)=x(r,fc^) + -^ 



7(3**^2 + k,k') 
327rG5 ■ 



(11.24) 



The counter term will cancel the UV divergences appear- 
ing in the expression of x and we will get a finite result at 
the boundary, i.e. limr->o X^'^'H^; ^m) ^il^ ^c finite. From 
the above analysis, we understand that one can get rid 
of the UV divergences appearing in the response func- 
tion by following the holographic renormalization tech- 
nique. But, an important observation is, this counter 
term does not add any finite contribution to the result 



it only cancels out the divergences. Thus, one can study 
the flow of the un-renormalized response function and ig- 
nore the divergences piece to get the finite contribution 
at the boundary. 



B. Second order transport coefficients from flow 
equation 

In this subsection, we compute the higher order trans- 
port coefficients by solving the flow equation (11.15 1 per- 



turbatively up to order lu^ and fc|. This is a non-linear 
flrst order differential equation. Now, the right hand side 
of this equation is proportional to w. Hence, to solve x to 
order uP , we can replace the leading order solution for x 



in the right hand side of equation (11.15 1. This simplifies 



the situation a lot as the non-linear equation becomes 
linear. Now, to leading order, x — ^^7 = — ig^g ■ There- 
fore up to order w^, we get. 



drXik^,r) = iuj 



9rr 

9tt 






T(r) 



+ 0(c.2,fcf). 



(H.25) 

The integration constant for the equation can be flxed 
form the boundary condition (11.18). Putting the value 



of the constant, the solution takes the form, 

1 



Jwx(/c;,,0) = lim 



r^o ge-TrGsr 



3q^r - 1) + Lj{3uj 



+r(a;(log(8)-3)+6i)) 
1 



167rGf 



+w" 



(l-ln2) 



1 
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O(^). 
r 



(11.26) 



Here we have chosen the four momentum in the following 
form 



{w,0,0,g} 



(11.27) 



This expression has divergence as r — > (UV diver- 
gence) and can be removed by adding suitable counter 
term (as explained in the last section). 

Comparing the finite piece of (11.26 1 at r ^- with 



the gen eric expansion of the retarded Green's function 
(|L2jQ we get. 



167rG5 
2 -In 2 
27rT ' 



ttT 



(11.28) 



The overall sign depends on the choice of coordinate. 



Here T = - . These results are in agreement with [TD] . In 
appendix [A| we briefly outHne the Kubo method to get 
this result. Thus, we see that, studying the flow equation 
of the response function we can compute the higher order 
transport coefficients perturbatively. Here, we present 
the results for the second order transport coefficients, 
but, in general it is possible to go beyond second order. 
At this point, it is not clear why only considering the 
boundary term from action (n.6 1 is enough to get the cor- 
rect results. In usual Kubo approach, one needs to take 
into account the Gibbons-Hawking term also. But as the 



tion. 



action (II. 6) has well defined variational principle, one 



does not need to add any Gibbons-Hawking term with it. 
In appendix [B] we show that the boundary terms coming 
from the original action and the corresponding Gibbons- 
Hawking action are exactly same as the boundary terms 
coming from the action (II.6) up to terms proportional 
to (fP' and pure divergence terms. The (j)^ terms do not 
contribute to any transport coefRcientfP] The divergent 
terms will get canceled by the proper counterterms and 
hence are not important for finding the transport coef- 
ficients. Thus it is clear that the effective action will 
give us the correct transport coefficients for the boundary 
plasma. This observation also holds for higher derivative 
gravity theorjj^ 



III. HIGHER DERIVATIVE CORRECTION TO 
FLOW EQUATION 

So far we have discussed the flow equation of two 
point correlation function of energy-momentum tensor of 
boundary theory whose gravity dual is given by Einstein- 
Hilbert action (two derivative action). But it is not ob- 
vious how to generalize this for higher derivative case. 
The proof given in ;27J was based on the canonical form 
( |II.6 | of graviton's action. In presence of arbitrary higher 
derivative terms in the bulk, the general action for the 
perturbation h^y does not have the above form (II.61. 
Rather it will have more than two derivative (with re- 
spect to r) terms like 4)'(t)'\ </>"^ etc. In presence of these 
terms it is not possible to bring this action into a canon- 
ical form (up to some total derivative terms). In this 
paper we consider generic higher derivatives terms in the 
bulk Lagrangian. We follow the prescription of [29] to 
construct an effective action "5'eff" for transverse gravi- 
ton in canonical form in presence of generic higher deriva- 
tive terms in the bulk. The effective action and original 
action give same equation of motion perturbatively in the 
coupling of the higher derivative terms. 

Let us consider a gravity set-up with n derivative ac- 



^^ They only contribute to pressure of the boundary theory. 

^^ In appendix IbI we have proved this explicitly for i?'"' gravity 



1 = 



1 



i? + 12 + a'i?(") 



(IH.l) 



where, TV-"^' is any n derivative Lagrangian. The metric 
in general is given by (assuming planar symmetry), 

ds^ = -{htir) + a' h^r\r))dt^ + '^"'^ , , 

Kir) + a' h^;'\r) 

+ -(l + a' /ii")(r))dx2 . (IH.2) 

Substituting the background metric with fluctuation 
(II. 3 1 in action (IH.l ) (we call it general action or original 
action) for the scalar field (p{r, k) we get, 



^=Y^ f^dr J2 A,Jr,k)cb^PHr,-k^(''Hr,k) 

(111.3) 

where, (j)^P>{r,k) denotes the p*'' derivative of the field 
(/>(r, k) with respect to r and p + q < n. The coefficients 
Ap^q{r, k) in general depends on the coupling constant a'. 
Ap^q with p + q > 3 are proportional to a' and vanishes 
in a' — )■ limit , since the terms (/>(^'(/)'^''' with p + q > 
3 appears as an effect of higher derivative terms in the 
action (IH.l). 



Up to some total derivative terms, the general action 
(II. 6 1 can also be written as. 



1 = 



1 
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leTrGs 



d'^k 

(2^ 



d*k 

J2ny 



^/2 



dr J2 Apir, k)(j)''P^ (r, -/fc)^^^) (r, k) 



p=0 



(for n even) 



dr J2 Ap{r, k)(t>^P\r, -k)4P\r, k) 



p=0 



(for n odd) . (III.4) 



However this action does not have canonical form. We 
write an effective action for transverse graviton in canon- 
ical form, 



ScS — 



1 



d'^k 



dr 
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+Anr,k)^{r,k)^{r,-k) 



A^''ir,k)<j>'{r,k)<t>'{r,-k) 



(IH.5) 



theory. 



with some unknown function A^^ and Aq^. We fix 
these functions by demanding that the equations of mo- 
tion obtained from the effective action and the original 
action are same perturbatively in a'. 

The generalized canonical momentum and equation of 
motion are given by, 

n™(r,fc)^2yl™(r,fc)0'(r,A:) 
(H™(r, fc))' = 2A^^{r, k)(b{r, k). (111.6) 



Once we find the effective action for tlie graviton, we 
follow the procedure in the previous section to obtain 
the flow equation for the boundary Green's function in 
generic higher derivative gravity. 

The boundary Green's function is given by, 



G^"{k^) = lim 



<Po{k)<j>o{-k) 



(111.7) 



which can be written using the definition of canonical 
momentum as, 



cr (fc^) = i™ 



n™(r,fc. 



r^o (j){r, kfj_) 



(III:, 



Let us define a response function of the boundary theory 
in higher derivative theory as. 



Therefore the flow equation is given by. 



drX ik^,r)^iuj 



9rr 

gtt 



r 

^^^(r,k) 



'-HD, 



(IIL9) 



r,k) 



where we define 



^HD 



(r,fc) = -2^r(r,fc^) 



T"^(r,fc) = 2X^(r,fc^) 



9rr 

9tt 

litt 

Qrr 
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(iii.ii) 

(III.12) 



This is the flow equation for two point correlation func- 
tion of energy-momentum tensor in presence of generic 
higher derivative term in the bulk action. Therefore 
integrating this equation from horizon to asymptotic 
boundary one can find the higher derivative correc- 
tion to the transport coefficients at any order in fre- 
quency/momentum. 

Like two derivative case here also we need to provide a 
boundary condition to solve this equation. The response 
function x^'^(^aii^) should be well-defined at horizon. 
This implies, 



r(fc.,r) 



EHD(r)THD(^) 



(III.13) 



here the horizon is located ai r — r^. 

One important point to mention here is that unlike 
two derivative gravity where x(A:^,?'/i) was independent 
of kfj^, x^^i^fj., Tfi) can in general depend on fc^. We will 
see this explicitly in the next section. Therefore the full 
momentum response at the horizon may not be able to 
correspond only to the zero momentum limit of boundary 
response in higher derivative theory. 

Like two derivative case, the response function in 
higher-derivative gravity theory also contains UV diver- 
gences. We need to add proper counter term following the 



holographic renormalization procedure to cancel these di- 
vergences. A little more thinking also says that in pres- 
ence of any higher-derivative term in the action the struc- 
ture of the counterterm remains same as (IL21I. Only 



the overall normalization constant depends on higher- 
derivative coupling. Thus, similar to the leading gravity, 
the counterterm in higher derivative gravity also cancels 
out the divergence and does not add any finite contribu- 
tion to the boundary response function. One can study 
the flow equation of the un-renormalized response func- 
tion and read off the transport coefficients from its finite 
piece. 



IV. EXAMPLES: STRING THEORY 
CORRECTIONS TO FLOW EQUATION 

String theory predicts next to leading order corrections 
to Einstcin-Hilbert action. These corrections are rele- 
vant at a distance comparable with typical length scale 
of the theory Is = \fa' . The short distance corrections 
to this action is described by supplementing this action 
by higher curvature terms. However, here we treat the 
stringy effects perturbatively i.e the coupling of higher 
derivative terms to be small. 

From the point of view of AdS/CFT the small a' cor- 
rection in supergravity corresponds to ^ correction in 
strongly coupled gauge theory in planar limit where A 
is the 't Hooft coupling. The precise dictionary between 
string length and 't Hooft coupling is. 



n 
a = 



4^ 



L AdS radius. 



(IV.l) 



In this section, we will consider two examples of higher- 
derivative terms coming from string theory and study 
their effects on flow equations and second order transport 
coefficients. 



A. Weyl^ term 



We consider the well known Weyl^ term. This term 
appears in type II string theory. Adding this term in the 
bulk action corresponds to -r372 correction in dual large 
N theory. The string theory correction to second order 
transport coefficients have already been computed in |32| 
using usual Kubo formula. Here, we show that one can 
obtain the correct result by studying the first order flow 
equation in higher derivative gravity without solving any 
second order differential equation for graviton. 

The five dimensional bulk action is given by, 

S^j^Jd'xV^(^R+12 + jW^^^) (IV.2) 



where, the coupling constant 7 is given by, 

1 



7 



;C(3)a' 



/3 



(IV.3) 



and 



From the regularity of x (^a" ^) ^^ horizon we get, 






-yhk' 



, 2G''^™"C,,™„C,'-^^C',,, (IV.4) 
and the Weyl tensors Cabcd are given by, 

Cabcd = Rabcd + T^igadRcb + QbcRad — QacRdb — gbdRca) 



12 



{gacQbd — gadgcb)R 



(IV.5) 



The background metric is given by [351 [37] VH 

- 15r {3r^ - Sr" + 5) jr^ ] dt^ 



,.._,(! -")•■.' 




1 (-285r'^ + 75r2 + 75)7' 



47.2 _ 4^4 



4- 4r2 



The temperature of this black hole is given by, 
T - ^ (1 + 157) . 



dr^ 



(IV.6) 



(IV. 7) 



and the horizon is located at r = 1. 

The effective action for transverse graviton has been 
computed in |29j . Here we will present the result only. 
The effective action is given by. 



I 



'off 



d-^k 



rdr 



AY {r,k)cl,'{r,k)cb'{r,-k) 
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+A^\r,k)<P{r,k)4,{r,~k) 



where, A^ and A^ are given in appendix O 
Therefore the flow equation is given by. 



(IV.. 



drX^ {kf,,r) =iuj, - 



gtt 



Si^'(r,fc) 



UJ^ 



where. 



(IV.9) 



S'^ {r,k) = ~2AY {r,k,) 



T^ (r,fc) = 2^S^ (r,fc^) 



9rr 



gtt 



gr 



(IV.IO) 



(IV.II) 



The explicit expressions for E and T can be obtained by 
using ytf" andy^S;^". 



^^ In this particular example we keep the extremality parameter ro 
explicitly. It would help us to write the relation between different 
transport coefficients. One can set ro to be one by time re-scaling 
what we have done in the next example. 



X^ (fcp,i) 



j:w'^(i)TW*(i) 



IGTrGs 47rGs 



(45r2 



llg^) 



(IV.12) 



Here, we see that unlike the two-derivative gravity, the 
horizon value of the response function depends on spatial 
momenta q (see Fig. pi). With this boundary condition 
we solve the flow equation up to order w^ and^^f^ (ignoring 
0{iijq^) term). Here we write the final resul 



ir I 
10 



itux"^ ik^,0) = -iil + 180j) 



16nG; 



-OJ 



(l-log(2)) 



-7(199- 66 log(2)) 
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+0{qu\ujq\q\u:^) . 



(IV.13) 



Comparing this result with (1.2) we get 



7r3T3 

K 



= 1 + 1357 + 0(7^ 



T^rT = 



ttT 

2 - log(2) 3757 



(1 - 1457) + 0(7^) 

+ 0(7^) . (IV.14) 



2tt 



47r 



These results are in agreement with 32 . It provides a 
non-trivial check to this approach of obtaining higher or- 
der transport coefficients from the flow equation (11.15 1. 



B. Four derivative term 

In this section, we will concentrate on the generic four 
derivative corrections to Einstein-Hilbert action. These 
terms arise in the effective action for the heterotic string 
theory. In fact, the complete super-symmetrized R^ cor- 
rection to effective Heterotic string theory is known and 
one way to obtain it is the super-symmetrization of the 
Lorentz Chern-Simons terms [351 ISS]- This terms also 
arises in the context of Type IIB string theory [101 [U] , 
where the theory is on AdS^ x X^, the compact space 
X^ being S'V^2- The dual theory is TV = 2S'p(7V)gauge 
theory with 4 fundamental and 1 antisymmetric traceless 
hyper-multiplets. This super-conformal theories arises in 
the context of iVD3-branes sitting inside 8 D7-branes co- 
incident on an orientifold 7-plane. In this case, generic 



{oj, 0, 0, q} and we ignore the UV divergence piece. 



four derivative B? correction comes form the DBI action 
of the branes. 

Here we compute the generic four derivative correction 
to the second order transport coefficients, the relaxation 
time T^ and k. We can choose the coefficients of the 
higher derivative terms to be the four-dimensional Euler 
density and get pure Gauss-Bonnet correction to these 
coefficients. 

The action 



1 = 



1 



IGTrGn 



x^/-g 



R+n + a'iPiR 



+/32i?M.p.i?^"''" + /33i?M-i?^" 



(IV.15) 



In particular for Gauss-Bonnet correction, /3i = 
l,/32 = l,/33 = —4. One can get rid of the Ricd? and 
Scalar'^ terms by a field redefinition and therefore all 
physical quantities should depend on the coefficient (^2 
only. Here, we prefer to work with the generic case as it 
would be easier for us the get the results for pure Gauss- 
Bonnet combination at every step. 

The background solution is given by |30| . 



ds^ = f{r)dt^ -f 



'^^ydr^ + Ux^ 



Zj.y3 J* 

where /(r) and g{r) are given by, 

/(r) = r - ^ - 2r(r^ - 
r ^ 

and 

2r(l0/3i-H(l-3r2)/32 + 2/33)a' 



(IV.16) 



2a 



(IV.17) 



9{r) 



l-r2 



3(r2 - 1) 



(IV.18) 



This is the background metric corrected up to order a' . 
We have fixed the integration constant such that the 
boundary metric is Minkowskian and the horizon is lo- 
cated at r = 1. The temperature of the black brane is 
given by. 



T ■ 



1 

■K 



10/?i - 5/32 + 2/33 



Btt 



-a . 



(IV.19) 



Similar to the Weyl'^ case, we can write the following 
effective action for this model. 



S, 



I 



eff 



d'^k 
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dr 



A^''{r,k)<t>'{r,k)<P'{r,-k) 



+A's''{r,k)^{r,k)^{r,~k) 



(IV.20) 



where, A^^ and A^^ are given in appendixpl Now, it is 
straightforward to write the corresponding flow equation 
(11.15 1 in this case. 



a,x^"(fc^,r)=zw./- 



9u 



X^Bfa,r)2 TGB(r,fc) 



(IV.21) 



where we define 

TGB(^,A;)-2^G^(r,fc^) 



gtt 



9tt 

9rr 



(IV.22) 



(IV.23) 



Now, the boundary condition (H.18) takes the following 
form, 



X^'^ik,,!) 



1 



[l + ((g2_8)/33-40/3i)a'] 



levrGs 

(IV.24) 
As mentioned earlier, we see that even in this case, 
the boundary condition depends on spatial momenta q 
through the coefficient P^. With this boundary condi- 
tion, one can solve the fiow equation (IV.21) and the 
solution is given by. 



iujx''''{k^,0) 



1 



levrGs 



i(l-(40/3i-t-8^3)a')w 



+- 



CO 



(l-log2)-f -^(130/3i(log2-l) 



-/32(51og2-2)-H26/33(log2-l)) 
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2 



l--(130/3i + 25/32 + 26/33)a' 



+0{quj',uq',q-',UJ-') 



(IV.25) 



From this expression we get the following transport co- 
efiicients. 



V 



levrGs 



(l-8(5/3i+/33)a') + 0(«'^)- (IV.26) 



This matches with results in gD) EH |13]- The higher 
order coefficients are, 

K = 4^ (1 - 10/32^') + 0{a'^) 
r^T = ^^ - ^^a' + 0{a'^). (IV.27) 



27r 



27r 



As we can see, the physical quantities r]/ s, k, t^^T only 
depend on the coefficient /32. In particular to Gauss- 
Bonnet combination, the corrections are, 

K= ^(l-10a') + C(a'^) 

ttT 

r^T = ^-^ ^a' + 0{a'^). (IV.28) 



(IV.27) and (IV.28) are new results of this paper. 



1. Exact result for Gauss-Bonnet black hole 

As we have done the above computation perturba- 
tively, the above expressions are valid only at order a'. 
But, one can consider the Gauss-Bonnet term exactly in 



10 



coupling. For pure Gauss-Bonnet combination the equa- 
tions of motion remain second order differential equation 
and hence it is easy to solve exactly to find the back- 
ground space-time. We solve the flow equation in this 
background exactly in coupling constant, and find the 
exact expressions for relaxation time r^ and k. In this 
section we briefly outline the result. 

The action and the solution is given by, 



1 



'-GB 



IGTrGg 



(fxy/^ 



R+12 



ds' = r'i-^^'^ 



dt^ + dx'^ 1 + 



'f{r) 



where. 



and 



2A„ 



/oo = lim /(r) = 



4A 



gb 



r4 



i-yr^4A 



gb 



2A 



■gb 



(IV.29) 



(IV.30) 



(IV.31) 



In this coordinate the boundary metric is 7]. We also 
consider only the '— ' branch of f± which corresponds to 
a non-singular black hole solution with non-degenerate 
horizon. 

The black hole temperature is given by. 



r = 






(IV.32) 



With exact GB, the effective action for fluctuation has 
a canonical form. Therefore we derive the flow equation 
for the response function (as we did in section pi) and 
solving this equation we get. 



T] 



1 



167rG! 



■(1-4A,,) 

2Agft (SAgb — 1) 



T^T = 



(1 - v/1 - 4A,b) (4Agfc 
1 



1) 



- 8X% + 12^/1 -4A<,bA 



47r(-l- 

+10Xgb - 2^1 - 4A<,fc - 41og(2)Agfc 

+ (1 - 4Ag6) log (-4Agfc + ^ 

-h(4A<,fc-l)log(l-4A36)-2 



gb 



4A 



gb 



log(2) 



(IV.33) 

One can easily check that up to first order in A^b, the 
results in (IV.33) reduces to the one in (IV. 28). In [44] 



the authors obtained the relation between second order 
transport coefficients and A^b numerically, however we 
are able to present the result exactly. 



As we have mentioned in introduction that the flow 
equation is a first order non-linear differential equation 
but one can reduce this equation to a second order linear 
differential equation. This second order differential equa- 
tion is related to the equation of motion for transverse 
graviton (gauge invariant excitations). Therefore we can 
use this equation to study causality violation in Gauss- 
Bonnet gravity. In [331 13S] it was found that to preserve 
causality of a conformal fluid there exists a bound on 
second order transport coefficients. 



r^T-2- >0 

s 



(IV.34) 



In Fig l3| we plot Tt^T —22 for our result and find the 
following bound on A^b which is in agreement with |44j . 



0.711 < A„b < 0.113 . 



(IV.35) 
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FIG. 3. Bound on A 
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FLOW EQUATION FOR CHARGED BLACK 
HOLES 



Electrically charged black holes in five dimensions have 
drawn a lot of interest in the context of AdS/CFT. The 
electric charge of these black holes are mapped to the 
global R-charge of the dual field theory. Because of the 
presence of the electric charges, the thermodynamics and 
the phase structure of these black holes are rather compli- 
cated and also interesting at the same time. There have 
been a lot of study of thermodynamics and phase tran- 
sitions of these charged black hole with different horizon 
topologies (see 0^ and references therein). 

The goal of the present section is to apply the 
AdS/CFT correspondence to understand how non- 
vanishing chemical potentials effect the hydrodynamic 
behavior of strongly coupled gauge theories. We study 
the second order hydrodynamics in two cases: (a) Generic 
R-charge black holes and (b) Charge black holes in 
higher-derivative gravity. 



11 



A. R-charged black holes 

We consider a conformal field theory witli conserved 
charge (density) in addition to energy and momentum. 
This is especially an interesting extension of the hydro- 
dynamics of the uncharged fluids. 

The second order hydrodynamics of charged fluid 
has been studied in [371 |3H]- They consider Reissner- 
Nordstrom black hole in five dimensions and found the 
effect of chemical potential on second order transport co- 
efficients in some limits of chemical potential. However, 
we consider generic -R-charged black holes with three 
(unequal) charges (chemical potentials) and find the ex- 
act expressions for second order transport coefficients in 
presence of three chemical potentials. As we have men- 
tioned in the introduction that solving the flow equation 
(of retarded Green's function of energy momentum ten- 
sor) we can only find two second order transport coeffi- 
cients whereas in [47l |48] all other second order transport 
coefficients have been reported. 

We consider R-charged black holes in five dimensions. 
A consistent truncation oi M = 8, D = 5 gauged super- 
gravity with 5*0(6) Yang-Mills gauge group, which can 
be obtained by S^ reduction of type IIB supergravity, 
gives rise to J\f — 2, D = 5 gauge supergravity with U{1)^ 
gauge group. The same theory can also be obtained by 
compactifying eleven dimensional supergravity, low en- 
ergy theory of M theory, on a Calabi-Yau three folds. 
The bosonic part of the action oi Af = 2, D = 5 gauged 
supergravity is given by [46]. We follow the notation of 



1 



levrGs 



d^Xy/^ 



R + V{X) 



-Gij{X)Fl^F>^''J - Gij{X)d^,Xid^'Xj 



+ : 



d^x e^'^P'^^ A^F.pF^^ (V.l) 
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where, AT^'s are three real scalar fields, subject to the 
constraint X^X'^X^ = 1. F^'s, which are field strengths 
of three Abehan gauge fields (I,J=1,2,3), and the scalar 
potential V{X) is given by, 

Fl, = 2d^,Ai^ 

G,,= 2'iiag[(Xl)-^(x2)-^(X3)-2] 



vw-^Y^i 



(V.2) 



The three-charge non-extremal STU solution is speci- 
fied by the following background values of the metric 



ds^ = -V.-'^'^ fk dr" + n^/^ {fk^dr^ + r^dQ 



3,fc) ' 



fk^k 



■nik 



r^n, 



H, 



% — H1H2H3 



(V.3) 



(V.4) 



as well as the scalar and the gauge fields 



X' 



H, 



V Qi 



The parameter k determines the spatial curvature of 
dilg f,: k ~ 1 corresponds to the metric on the three- 
sphere of unit radius, fc = - to the metric on R^. Hy- 
drodynamic approximation is valid only in the case of a 
translational-invariant horizon, in our case we set fc = 
and 

d^l.Q -^ {dx^ + dy^ + dz^) ■ 

Replacing the radial coordinate r -^ J'o/nA'i where rg is 
the largest root of the equation /(r) — 0, the background 
solution in this new coordinate is given by, 

dsl^ -n-^'^^^^ I dt^ +'H^/^-^dr^ 
r 4/r^ 



n 



m i^nr 



[dx^ + dy"^ + dz^) 



(V.6) 



/(r)=H(r)-r2]J(l + «,), H, = 1 + k, 



(V.7) 



where n'^s are chemical potentials and 

To - ro/^ . (V.8) 

The scalar fields and the gauge fields are given by 



A' = 









where, 



AZi JJ(1 + Ki) 



1/2 



(V.9) 



(V.IO) 



The Hawking temperature of the background (V.6 1 is 
given by 



„ _ 2 + Kl + K2 + K3 - K1K2K3 
"^ 2^/(1 + Kl)(l + K2)(l + K3) 



(V.U) 



We perturb the xy component of background metric 
and the action for transverse graviton is given by. 



S'cff — 



1 



d'^k 



dr 
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+A^{r,k)ci,{r,k)ci,{r,-k) 



A'^{r,k)(j)'{r,k)(t)'{r,-k) 



where, 



A? 



rUir) 



(V.12) 



(V.13) 



12 



and 



-^0 



Q _ ^0 
4^2 



H1H2H3 

~1W~ 



q 



(V.14) 



Therefore the flow equation is given by, 



drX^ikf^r) = iujJ- 



Qrr 

9tt 



i:^{r,k) 



(V.15) 
Solving this equation perturbatively in uj and q we get 



x'^iK^r) 



ro'n.(i 



U/2 



leTrGs 
'2uj 167rG5 



^aa+^o 



16^G5V4^K + (l + 5« 

^^^ 1 + ^« - 2P, - V4P« + (1 + ^«)2 



hi 



i + s^- 


- 2rS, - V4P. + (1 + ^.)2 


i + s^- 


- 2rP, + V4P. + (1 + 5«)2 


i + s^ 


- 2S^ + v/4P« + (1 + S^)^ 



+0{qLJ^Ujq^,q^Lu'') 



(V.16) 



where, 



P. 






H' 



(V.17) 



Computing the response function at the boundary 
(throwing away the divergent piece) we get the foUow- 
ing transport coefiicients. 



V 
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Y[a + ^^) 



1/2 



(V.18) 



and 



T-irT ^ 



77 1 + SJ2-PJ2 
2 + S,-P,. 2_ 



'3 + 5„ 



n.(i 



v/4P« + (1 + 5« 



In 



V^ 



+ (1 + ^«) 



3 + 5„-v/4Pk + (1 + 5« 



(V.19) 



These are the new results in this paper. It is easy to check 
that for Ki — >■ limit we recover the results in section [ll] 
To complete the discussion on the second order trans- 
port coefiicients for _R-charged black holes one should find 
the flow of Green's functions for two point correlation 
functions of _R-currents. As we will mention in section 



VI 



that in presence of finite charges (or chemical poten- 
tials) it is very hard to solve the Riccati equation even 
perturbatively in lu and q. We find it very difficult to 
get any analytic solution for i?-current Green's function. 
However we consider a simple model in section |VI| and 
study the flow of i?-current Green's function numerically. 



B. Charged black holes in higher derivative gravity 

In this section, we will study five-dimensional gravity 
in presence of a negative cosmological constant and cou- 
pled to C/(l) gauge field. The model has been studied in 
[551 [Sni [5T], the action is given as. 



^=16^/^'"^^ 



1 9 
i? + 12- -F^ 



C.a 



bcdt 



^aFbcFde + «' f^-lRabcdR^' 



abed 



+C2RabcdF"-'r^-' 

+C5e 



bcde AT} p. 






CiF^ 



(V.20) 



Here, F^ = FabF^'^F^ = FabF^" F^dF'^'' , and the AdS 
radius is set to unity. The action includes the Chern- 
Simon term and also a generic set of four derivative 
terms. All the four derivative terms will be treated 
perturbatively in our computation and here a' << 1 
is the perturbation parameter. In |33j . it was shown 
that, within perturbative approach, after using field- 
redefinition, this is the most generic four derivative action 
that one can write down. In this section, we will closely 
follow their work. The background metric and the gauge 
field in presence of these higher derivative terms have the 
following form, 

ds^ = ~r^f(r)dt^ + -^^dr^ + r^idx^ + dy^ + dz^), 



A == h{r)dt, 



(V.21) 



where, 



/(r) = /o(r)(l + a'F(r)), 
<7(r) = /o(r)(l + a'(F(r) 
h{r) = hQ{r)+a'H{r). 



G{r))), 



(V.22) 



Here fo{r),goir) and ho{r) are the solution of the back- 
ground in absence of the higher-derivative terms in the 
action and they are given as. 



/o =30 



1 



1 



Q' 



ha = VSQ [ ^ ^ 



(V.23) 



Here, Q is related to the physical charge of the system 
and tq is the position of the horizon. From (V.22), it is 



clear that even in presence of the higher-derivative terms, 
the horizon remains at tq. The higher-derivative cor- 
rections to this background are given by the functions 
F{r),G{r) and H{r). The form of these functions are 
given in |33j . We would not write those expressions and 
refer the reader to that paper. 

Using the flow equation, we will study the higher order 
transport coefficient of the plasma theory dual to this 
gravity model. For this, we will write the effective action 
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for the metric fluctuation in (II.3), as we have done in 
previous sections, 



•^cff — 



1 



d^k 



dr 



A?''{r,k)<f>'{r,k)c^'ir,~k) 
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+^^^(r,A;)0(r,A;)0(r,-fc) 



(V.24) 



where, Ai^ and Aq^ are given in appendix^ The cor- 
responding the flow equation (11.15) for this case with 
the coefficients Ai^ and Aq^ is, 



drx''''{k^.,r)=iLuJ- 



9tt 






(V.25) 



where we define 



S^"(r,fc) = -2^r(^,^M) 



9rr 

9tt 



CB/ 



T^"(r,fc) = 2A^"(r,fc^). - 



9tt 



9r 



(V.26) 



We solve this fiow equation to find the effect of higher 
derivative terms and chemical potential (or charge) on 
transport coefficient. Here, we present our result for 
small Q only, though it is possible to find the results 
for any Q. 



The boundary condition (11.181 will take the following 
forrrp^ 



Xikf^^ro) 



IQ-KGr, 



1-a 



,24ciQ2 



(V.27) 



In this case, the horizon value of the response function is 
independent of momenta. With this boundary condition, 
we can solve the fiow equation (V.25), and the solution 
is given by. 



It is easy to read off the transport coefficients from this 
expressions. 



'0 



167rG5 

■n 



1-a' 



,24ciQ2 



0{Q') 



t:T 



1- 



2rl 



a'{ 10ci-^(37ci-48c2; 



-O(Q') 



T^T 



2 -In 2 Q^{5-31ti2) 



2n 



Anr^ 



^^""^ +^(-16c2 + 5ci(ll-41n2)) 



2tt iirr^ 



(V.29) 



where, the temperature T of the system is given by, 



T = 



ro 



1- 



24 



5ci^ 



g(31ci+48c2) 



-0{Q'). 



3 V"^ ' 2r6 

(V.30) 

We see that the first order as well as the second order 
transport coefficients coming from retarded Green's func- 
tion of energy momentum tensoiP^ only depends on two 
coefficients Ci , C2 . This feature was observed in 33] for 
entropy density s and first-order transport coefficients ry. 
The coefficients C3,C4, which parameterize couplings in 
the four point function of the dual U{1) current does not 
play any role in these hydrodynamic coefficients. They 
should be important for the computation of conductivity, 
which comes from the Green's function of the boundary 
_R-current. Two other coefficients (, C5 also do not ap- 
pear in the expressions. One can find a magnetic brane 
solution of the action (V.20) like [S^- In that case it 



would be interesting to find the effect of magnetic field 
on transport coefficients. 



iwx(fc^,oo) ^iuj 



'0 



167rG5 
2 



1-a 



,24ciQ2 



— w 



'0 



327rG5 



(l-ln2) 



Q' 



(3 -In 2) 



--(ci(2-51n2) 



-48c2(2-ln2)) 



2rl 



(ci(35-581n2) 



grVg 



1-- 25ci + ^(32ci + 24c2) 



327rG5 



^'^ Note that we are working in a different coordinate where r — > 00 
is the boundary, therefore we choose the positive branch of the 



boundary condition dll.lSll 



VI. FLOW OF RETARDED GREEN'S 
FUNCTION OF BOUNDARY R CURRENT 

Finally, in this section, we study the fiow of retarded 
Green's function of boundary i?-current, 



G^jik) = -i / dtd^xe'''-mMx), Jj{0)]) {VIA) 

where Jf^(x) is the CFT current dual to a bulk gauge 
field A„. 



(V.28) 



^* It would be interesting to study the flow of retarded Green's 
function for boundary i?-current. We found it to be difficult to 
get any analytical solution for response function in presence of 
finite chemical potential and higher derivative terms. However, it 
would be nice to know the higher derivative corrections to other 
second order transport coefficients appear in R current 1471 148| . 
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In hydrodynamic approximation one can express the 
current in powers of boundary derivatives. Up to first 
order in derivative expansion it has the following form, 



Ai 



J, = -kP^d„^+m, + 0{d') 



(VI.2) 



where, k and 57 are two first order transport coefficients, 
/i is chemical potential, T is temperature and 






(VI.3) 



The expression of J^ up to second order in derivative 
expansion can be found in |47l 148] . From conformal in- 
variance of the theory it is possible to write all possible 
second order transport coefficients appear in the expres- 
sion of J^. However, like energy momentum tensor, from 
the expression of retarded Green's function it is not pos- 
sible to compute all the transport coefhcients that appear 
in different order of derivative expansion. 

In this section we study the flow equation of retarded 
Green's function of boundary i?-current. Unfortunately 
we find it difficult to solve the flow equation analyti- 
cally to extract any transport coefficient. We present our 
calculation how to write the flow equation for retarded 
Green's function of i?-currents in presence of generic 
higher derivative terms in bulk Lagrangian and some nu- 
merical results. 

We start with Einstein-Maxwell action 



^=i6k/^'-^(^+^^4^^ 



temperature of the black hole is given by. 



Solution is given by equation (|V.21|) with a' = 0. The 

(VI.5) 






TT 



and the chemical potential is given by, 

1^= —J2-- 



(VI.6) 



For technical advantage we write the metric and gauge 
fleld in a different coordinate. We change the radial coor- 
dinate r — >■ -^. In this coordinate the metric and gauge 
field is given by. 



ds' 



rlU{r),. 



-dt^ 



4r2[/(r) ^ r^ ' 



Atir) = E{r) 



(VI. 7) 



where. 



U{r) = (l-r)(l + r- 



Q\^ 



(VI.8) 



We turn on small fluctuations for x component of gauge 
fields. Since the At component of the bulk vector is non- 
vanishing in this background, the perturbations A^ can 
couple to the tx component of graviton Therefore we also 
need to consider small metric fluctuations for components 
gtx- Writing them in momentum space 



Ax{r,x) 



d^k 

d^k 

(2^ 



e'''-'''Ai{r,k) 



(VI.9) 



However, there exists a constraint relation between A^ 
and gtx ■ We use this relation to replace gtx from equation 
of motion of Ax ■ 

The on-shell action for gauge field fluctuations are 
given by. 



Sa^ 



1 



d'^k 



167rG5 J (27r)4 
( ' 



-TlU{r)A'^{r,k) 



-/It, 



(VI.IO) 



\ArU{r) Ar ^ 
The current corresponding to Ax fluctuation is given by, 

5Sa 



Jx{r,k) 



6A'^{r,k) 

-2rlUir)A',{r,k). (VI.ll) 



(VI. 4) 'pijg equation of motion for Ax{r, k) is given by. 



{U{r)^,{rM' = 



1 

47-0 \rU{r) r 
-E'{r)(j)'{r,k) 

1 / u;^ g' 

Atq \rU{r) r 



Ax{r,k) 



Axir,k) 



+ 



rE'{rf 



Ax{r,k). 



(VI.12) 



Using the constraint relation (coming from rx component 
of Einstein equations) , 



rE'{r)Ax{r,k) 

<P (r,k) = 2 ■ 

'0 



and the equations of motion one gets. 



(VI.13) 



J'xir,k) 



UJ 



2 \rU{r) r 
^2rE'{rfAx{r,k). 



Ax{r,k) 



Next we deflne a response function 

Jx{r,k) 



a{r, k) 



iujAx{r,k)' 



(VI. 14) 



(VI.15) 
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Taking the derivative with respect to r and using the 
equation of motion we find the flow is given by, 



a'{r,k) = 



2roU{r 



a(r,kf 

q^U{r)\ ArQrE'{rfU{r) 



1- 

r \ uj 



(VI.16) 



From the regularity of the response function at the hori- 
zon we find the boundary condition is given by, 



cr(l,/c)2 = l. 



(VI.17) 



With this boundary condition one can integrate this 
nonhnear equation to find finite frequency response of 
boundary Green's function. 

In presence of generic higher derivative terms in the 
Lagrangian the on-shell action for fluctuation Ax may not 
have canonical form like (VI. 10). In that case one has to 



write an effective action for A^ like transverse graviton 



The effective will have the same form as (VI. 10) only the 



coefficients will depend on coupling constant of higher 
derivative terms. 

We conclude this section by presenting some numeri- 
cal solutions of the ffow equation (VI.16) for both non- 
extremal extremal black holes in fig. |4j 




FIG. 4. Flow of a for non-extremal and extremal black holes. 
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Appendix A: Second order transport coefficients 
from Kubo formula 



d^xy/^[R + 12] 



(A.l) 



now consider the action with solution given in section ITTj 

For well defined variation of this action we need to add 
a Gibbons-Hawking boundary term. Also, requiring the 
on-shell action being finite at boundary, we have to add 
counter-terms following usual approach of holographic 
renormalization. They are as follows: 

Sgh = „ „ / d'^'xy/^JC, 



ScT = 



SttGs 
1 



1 



16.G../^'^^[^+2^] (^-2) 



where 7 and TZ are boundary metric and Ricci scalar 
(constructed out of 7) respectively. 
We consider the following metric perturbation, 

9.y = 9^°J + hxyir, x) = g'i^ (1 + e<i>(r, x)), (A.3) 

where e is an order counting parameter. We are in- 
terested in quadratic on-shell action for this transverse 
graviton $(r, x). Let us define the Fourier transform, 

d'^x 



(j)(r, k) 



-ik.x 



{2^Y 



<^{r, x) 



(A.4) 



and k — {— w, fcj. S ubstituting this fluctuation in action 



(A.l), we get (B.l ■^ Now, we rewrite this action (A.l) 



as equation of motion piece (which vanishes on-shell) and 
boundary term[^ Thus, on-sheU S, Sgh, Sct become, 



S, 



OH 



S, 



16ttG5 

1 

SttGs 

1 






=5 (27r)4 
d'^k 



CT 



LGH{Hr,k)), 
LcTiHr,k)), (A.5) 



le^Gs Jr=5 (27r) 
where we define, 

L((/)(r, fc)) = a2(r)(/)'(r,fc)(A*(r, fc) 

a4(r) — a5'{r) , , , , ,^ , , , 

+ 2 '^^'^' ^^^''' ^' 

LGH{cl^{r, k)) = 2(5l)0(r, k)(t>* (r, k) + 2(52)^'(r, fc)^ (r, fc), 

Lcri^ir, k)) = (co -I- c^uj^ + C2q^)(j){r, k)cb* (r, k). (A.6) 



The coefficients 02,04,05 are given in (B.2) and other 
coefficients are. 



1 



51 



Co = -3 



32 



VT^^ 



ci 



4r 



Cl 



4r\/r 



(A.7) 



In this appendix we re-derive second order transport 
coefficients k and r^r from usual Kubo approach. Let us 



^^ as , ag are zero here 

^^ We ignore contribution from liorizon as [5] 
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The retarded Green's function is defined at asymptotic 
infinity as, 



G^{k) = 2 fim 



(L + LcH + LcT)\on-shell 

r^o (j)a{k)(l)o{-k) 



(A., 



Here , 4'o{k) is the boundary value of the fluctuation 
(A.4). The retarded Green's function G^ is a function 
of boundary momenta k^ — (w, 0,0,(7). Now the lead- 
ing action and the Gibbons-Hawking action get diver- 
gences from 4>'{r, k)(f)*{r, k) and (j){r)(l)*{r, k) parts. Both 
these divergences get canceled by the counter-term ac- 
tion which is always proportional to only (f)(r)(f)*(r,k). 
In this case of leading Einstein's gravity, it is even more 
simplified. Divergences coming from (f)(r)(f>* (r, k) piece 
of leading and Gibbons-Hawking action gets canceled by 
momentum independent piece of Gounter-term action. It 
turns out that there is a cancellation among the corre- 
sponding coefficients as, 



lini(-(a4(r) 



a5'(r)) + 2.gl 



1 



Co) = 



(A.9) 



i.e. the final contribution from - 



' (r, k) piece is only a 



finite number ^. As the graviton fluctuation (j){r, oj, k) = 

(j)o{l + F{r,uj,k)) and moreover liuir^o F{r,uj,k) — 0, 
we see that the (t){r)(j)* {r) term above only contribute to 
pressure (the lo independent piece of G^). It would never 
contribute to any transport coefficient. 

Also, the divergences coming from 4>'{r^ k)(t)*{r, k) piece 
of original and Gibbons-Hawking action, get canceled 
with the piece of the counter-term proportional to ur^ , q^ 
(ci and C2 are purely divergent at boundary). Here, 
the situation is more subtle, as there is no cancellation 
among the coefficients. One actually needs to put the 
solution of 0(7-, w, k) to see the cancellation. 

The overall lesson from this detailed analysis is that 
counter-term only cancel the UV divergences in usual 
holographic renormalization process and at most con- 
tribute to pressure of the boundary plasma. It has no 
effects on any transport coefficients. In [10 , the author 
have computed second order transport coefficients for the 
plasmas dual to leading Einstein's gravity following this 
usual approach. The results are as follows, 



_ 2 - In 2 _ 77 



(A.IO) 



These results match with the one we obtained in (11.28) 
by solving the fiow equations. 



Appendix B: Equivalence of Boundary Terms 

In this appendix we will show explicitly why the trans- 
port coefficients computed from the original action and 
the effective action are same, even for any higher deriva- 
tive theory. It was already noticed [H], that the two 
would give same first-order transport coefficient rj with 



a suitable choice of the overall normalization constant. 
Here, we show that, not just the first order transport co- 
efficients, rather any higher order transport coefficients 
computed from the original action and the effective ac- 
tion are same. 

We consider a general class of action for (f) which ap- 
pears when the higher derivative terms are made of dif- 
ferent contraction of Ricci tensor, Riemann tensor, Weyl 
tensor, Ricci scalar etc. or their different powers. Since, 
all these tensors involve two derivatives of metric they 
can only have terms like dadb^{r, x) and its lower deriva- 
tives. Therefor the most generic quadratic (in $(r, a;), 
in linear response theory) action for this kind of higher 
derivative gravity has the following form (in momentum 
spacerH 



S = 



1 



cf^k 



dr 



al{r)(j){rf + a2{r)(j)' {rf 



IGTrGs J {2tiY 

+aA{r)(t>{r)(t>' {r) + a' a6(r)0"(r)(/)'(r) 

+a ai{r)(j)" [r) + a^{r)(f){r)(j)" [r) 



(B.l) 



where. 



a\{r) 
a2(r) 
a4(r) 
a5(r) 



4/^3 — 4r^ 



«'/2(r) 



-3r + - -I- a' h2{r) 
r 

-^-2 + a' g2{r) 

4 
-4r H h a' j2{r) 



(B.2) 



and a3{r) , a6{r) , j2{r) , g2{r) , h2{r) and /2(r) depends on 
higher derivative terms in the action and hence are com- 
puted purely from the background solution with a' — >■ 0. 
Among these coefficients a3 is special, as, it couples to 
(j)"'^. All four derivatives act on the graviton fluctu ation 
and thus a3 only depends on metric functions in (II. 2) 
and there r-derivatives. It is easy to convince ourselves 
that a3 ex r{r'^ — l)^/(r, a'), where /(r, a') is a function 
that depend on the higher derivative terms and finite 
(constant or 0) at the boundary r — j' 0. Now let us write 
the effective Lagrangian as follows. 



Seff - 



1 + a'T f d^k 



167rG5 J (27r)4 
+a'(62(r)0(r)2- 



dr 



4r{r^-lf(f>'{ry 



4^2 (-^2 _ l-j 



bl{r)<l)'{ry 



(B.3) 



Demanding that the equation of motion (up to order a') 
of (f) derived from the original action and the above action 



In all the expressions we have omitted k dependence of < 
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are same we get, 



are given by, 



61(r) 



1 



2r(r2-l)- 



r[(-4r3-12r + a;2)a3(r) 



62(r) 



+ {r^ - l){2Kr^ - a6' (r)r^ - Akt^ + 2a3\r)r^ 
+2(r2 - l)h2(r)r - 2{r^ - l)j2(r)r + a6'(r)r 
+2K + 2a3'(r))] (B.4) 



(w'' + 144r'^w'')a3(r) 



16r2 (r2 - 1)' 
+4(r2 - 1) ( - 4r2/2(r)(r2 - 1)3 



+{{uj^K, - 2r2(r2 - l)j2" (r))(r2 - 1) 
+2r^g2'{r){r^ ~ 1)2 + ru^ a3" {r)) (r^ - 1) 



+(1 - Ilr2)a;2a3'(r) 



(B.5) 



The boundary terms coming from the original action (af- 
ter adding Gibbons- Hawking boundary terms) are given 
byR 



S' 



1 



+ 0(r)2+r<?!)'(r)0(r) 



167rG5 J (27r)4 
^'^^+a'(^lg2irmrr-lj2'irmr? 

+(/i2(r)-j2(r)-"^)(/)'(r)0(r) 

a3'(r) ((/)(r)cj2 + 4 (r" - l) cj,' (r)) <j){r) 
4^(j.2 _i)2 

a3(r)(6r<?!)(r)(/)'(r)a;2) 

4r(r2- 1)3 

a3(r) ((r2 - l) (Sr^ + 24r ~ lj^) (/)'(r)) 0(r) 

4r (r2 - 1)3 

a3(r)0'(r) (0(r)cj2 + 4 (r* - l) (/)'(r)) 

4r (r2 - if 



-a3(r)0'(r) 



2r{r^-iy r{r^-lY 



(B.6) 



And the boundary terms coming from the effective action 



s: 



seff 



1 r <Pk 

167rG5 j (2^ 

a' 

2r(r2-l)2 



r 1 (f){r)(j) [r) 



'/-^^^3 



(r)(2r(r2-l) +(-a6'(r)r 



+2a3'(r)r2 + 2 (r^ - l) h2{r)r - 2 (r^ - l) j2(r)r 
+a6'(r)r + 2a3'(r))(r2- l) 

+ (-4r3 - 12r + lo^) a3(r))0'(r) j 



(B.7) 



Now, it is interesting to compute the difference between 
these two boundary terms and the resuh i^\ 



nB qB _ -^ 



d^k 

(2^ 



0(r)2 



Hr? 



-j2'(r)0(r)2 



+a'(^-52(r)0(r)2 

a3'(r)a;2^(r)2 a3(r) (6rcj2) 0(r)^ 



4r (r2 - 1) 



4r (r2 - l)-" 



(B., 



The term proportional to a3 in the parenthesis of (B.8) 



vanishes at the boundary whereas the term proportional 
to a3' gives a pure UV divergent piece and a vanishing 
piece, due to the property of a3 mentioned above. This 
is true irrespective of the choice of the higher derivative 
terms. Thus, we see that the two boundary terms differ 
only by terms which are either purely divergent or of the 
form g(r)(/>2, where g{r) is any function of r. The diver- 
gent terms would get canceled once appropriate bound- 
ary ter ms are added (which has been discussed in sections 
O and III ) . The terms proportional to cfP can only con- 



tribute to pressure of the boundary theory and are not 
important for the computation of transport coefficients 
of the boundary plasma. Thus we see that, it is obvious 
that the transport coefficients coming from the original 
action and the boundary action are same. 

Here we have considered only i?*^") gravity theory. A 
more rigorous proof is required for theories involving co- 
variant derivatives of curvature tensors and scalars. 



^^ There was a sign error in |29| 



It has been shown in i29i that F = 0. 
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Appendix C: Functions appeared in String theory 
corrected action 



Appendix D: Functions appeared in four derivative 
action 



A^^{r,k) 



g2 ^^2 „ ;l^ ^_ ^2 
4^2 (^2 _ Ij 



12,2(,2_i)y9'(2/33(13r^-3r^^-13) 
+130(r2 - l)/3i + (-36r'' + 25r^ + ll)/32) 
+a;2((6r2 _ -^^^^^ _^ ^3Q^^ _^ 26/33) ) (D.l) 



and 



A^ (r,fc) 



rg(g2(r2 - 1)+^^) 
4^2 (y,2 _ 2) 

1 



4(,2_i)^7(49V((245r^ 
-407r2 + 162)ri^ + lOrw^) 
+3 (221r'^ - igir^ + 25)wV^) 



(C.l) 



^^«(r,/fc)=r 



(r2 - 1) ((I8r2 - 13) /32 + 110/3i + 22;93) a' 



3r 



(D.2) 



Appendix E: Functions appeared in 
Higher-derivative Charged black-hole action 



and 



< (r,fc) 



r'-l]r^n 



-3r(r2-l)7r^(43r'' 
-47r2 - 25)ro + 16gV^ 



(C.2) 
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CB 



cia'(llr8(a;2 _ ^2) _^ ^4^4(25^2 _^ g^2) _^ 36^2^8) 



Q^ 



2r2(r2-r2)(r2 + r2)2 



67-3^^4 _ ^4 -J 

3, ,2 j_ _^ 

3r5 



r"a;' + tttt I ci[2%q^rlr'^ 



- r^{-i&q^ + 127^2) - Arlr^{7q^ - 30;^) - 8q^ 

+ 4r^r''(llg2 + iu?)] - 2Ac2[rlr^ {2q^ - iu?) - 2q^rlr 



,2„6„2 



+ roV(2g2_3c.2)_2gV8+4r« 



0{Q^) (E.l) 



and 



,CB_1. 4 ^5^ cic.- (I3r^ - 18rg) (r^ - r^^) 



+ 



Q2(^2_^2) 



2rrg 



3r4 



24c2(4r''-3r^r^-3r^) 



+ci(101r* - 156r^r^ - 120rc;) 



0{Q^). 



(E.2) 
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